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Abstract 

We derive a simple rule to determine surface plasmon energies, based on the geometrical prop- 
erties of the surface of the metal. We apply this concept to obtain the surface plasmon energies in 
wedges, corners and conical tips. The results presented here provide simple and straightforward 
rules to design the energy of surface plasmons in severals situations of experimental interest such 
as in plasmon wave guiding and in tip-enhanced spectroscopies. 

PACS numbers: 
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I. INTRODUCTION 



Plasmonics has become an important branch of nanooptics that allows for the manipu- 
lation of light in the nanoscale thanks to localized surface collective modes so called surface 
plasmons. Surface-plasmon propagating in plasmonic circuits^ may act as interfaces in 
electrooptical devices, and localized surface plasmons play a key role as the mechanism of 
field-enhancement in field-enhanced spectroscopies assisted by metallic nanoantennas^. The 
capacity to tune the optical response in metallic systems results of crucial interest in many 
processes of physics^, chemistry^, biology^ and medicine^ where metallic systems act as 
active hosts responsible for the existence and efficiency of the processes. 

The energy of surface plasmon is known to depend on the metallic material, the geometry 
of the system, and the environment. Many different geometries have been studied from the 
simple metallic sphere or spheroids^ treated analytically in the beginning of last century to 
complicated coupled systems tackled with use of sophisticated numerical methods to solve 
Maxwell's Equations^!. Here we derive a simple rule to determine surface plasmon energies 
in direct connection with the metallic shape. Laplace equation is expressed in terms of a 
surface charge density at the boundaries of the metal surface, and solved using a standard 
method of integral equation. This method allows for straightforward connection between 
the surface plasmon energies and the solid angle at the metal-vacuum interface. 

The simplicity of the concept allows for direct estimation of plasmon energies in common 
situations of interest such as in a wedge waveguide with plasmons running along a metallic 
wedge, or at the tip of a metallic cone where surface plasmons are localised at the tip apex. 
This piece of information permits to estimate and design the energy of surface plasmon that 
can be later fine-tuned to engineer and control their functionality. 

II. EXTENDED THEORY 

Here we investigate the electromagnetic modes localized at the surface S separating two 
different media I-II, characterized by two dielectric constants ei and 62 respectively. The 
normal vector n s at the interface position s has been chosen arbitrarily to point towards 
medium II, see FIG. 1. In this way, we shall call "internal" the variable defined in the part 
of space corresponding to medium I, and "external" the variable defined in the part of space 
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Medium I] 




FIG. 1: Interface between two regions with different dielectric constants e\ and 62 bounded by a 
surface S. 

corresponding to medium II. For instance, t\ will be denoted in the following by e int and e 2 
by e ext . 

We want to obtain the surface plasmon modes of this system. These modes can be de- 
scribed by the Laplace equation for the quasi-static potential u(r) of the self field, everywhere 
except at the surface S where we have the following equation 



VVr) 



-o r 



(1) 



with er(r) the surface charge density, i.e. different from zero on the contour surface S. Let 
us note by x the position vectors located on the surface, therefore cr(x) is related to the 
normal derivatives of the potentials internal Ui n t and external u ex t to the surface S by the 
following relational 

, \ dui n t du ext 



dUa 



(2) 



The continuity equation for the self-fields at the boundary gives the following condition 
between the normal derivatives of the potential at any point of the surface S 



duj. 
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From the previous equations, we can recast Eq. (1) as follows 

«(x) = / Cr(x, x / )cr(x / )(iS'(x / ) 
Js 



,-, i Oil dll ex t i , 

- dUi n t 



= (1 - e int / eert ) / G^dS(x') (4) 

where G(x, x') = 1/4tt|x — x'| is the Green function solution of the homogenous Laplace 
equation 

V 2 G(x,x') = -<5(x-x') (5) 

where 5 is the Dirac delta function. 

Now considering a point x = x G on the boundary S with possible singularities, and 
applying eqs. (A5) and (A6) of the appendix, i.e. evaluating the integral in Eq. (4) in 
Cauchy Principal Value sense, we obtain the following integral equation^ 

^( x d) e ext 



U(X ] 



4tT 6 ext — £i n t 



u( 3 o(~(x o ,x0)d5(x / ) (6) 



where fi(x ) is the solid angle sustained by the surface at the point x . 



Let us measure solid angles in units of 47r(f2 = and define 5£l = £1 — 1/2 as the 
deviation from a smooth surface. Notice that SQ < for a surface with a protrusion 
(convex) and 5Q > for a concave surface. We can now write 

/ d£l(xQ, x')u(x') = Am(x ) (7) 
Js 

where we have used the fact that J s (— fn) dS(x!) = , and consequently rff2(x ,x') is 
the normalized elementary solid angle sustained by dS(x') viewed from x . To the best 
of our knowledge, the surface integral equation allowing to obtain the surface modes of 
the corresponding system was never written as a solid angle integral equation Eq.(7). This 
equation is pointing out the purely geometrical nature of these modes. We note that equation 
(7) is an exact expression, that complements recent numerical approximations^ for smooth 
surfaces, and puts them into a more general context that includes surfaces with singularities. 
Finally we interpret equation (7) as an eigenvalue equation where 

A = [AqM + SCl) (8) 



and we have denned the frequency factor 

AoM = J;" + _y , = 1/2 - ( W / Wt ) 2 (9) 

In the last line we have introduced the generic case of a Drude metal for tint with bulk plasma 
frequency u h and vacuum value for e ext . Notice again that the right hand side of Equation 
(6) is entirely a geometric object. Furthermore, since the integral is evaluated in the Cauchy 
Principal Value sense its main singularity is already present in the factor Sfl = Cl — 1/2. 
Now let us consider the A = eigenvalue 

e«rt(l + e int Q = (10) 
In the Drude limit of Eq.(9), we obtain the following for the plasmon energy uj 

(u/u b ) 2 = n. (ii) 

III. RESULTS AND DISCUSSION 

The derivation of the previous section has shown the crucial connection between the 
shape and the spectral properties of the surface plasmon modes. In this section, for some 
selected geometries, we demonstrate the existence of the A = eigenmode. Afterward, we 
apply Eq. (10) corresponding to the A = case for metallic-vacuum interface, using a Drude 
dielectric response for the metal, and obtain the energies of these modes. 

A. A = mode for a semi-infinite metal 

The planar interface is a well known geometry for its analytical solutions of the Laplace 
equation. For a semi-infinite metal, with xo on the planar surface, any elementary solid angle 
of the surface, view from x vanishes (dii(x ,x') = 0). Consequently any function defined 
only in the interface, will be solution of the integral equation (7), with eigenvalue A = 0. 
For a planar surface, as we are dealing with a smooth surface, we have Cl = 1/2. Therefore 
applying Eq. (11), we obtain the well known infinitely degenerated surface plasmon mode 
of frequency oo s = ujb/VZ- 
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B. Modes localized in one dimensional surface singularities 

Now let us imagine a transformation, which generates from a smooth interface an interface 
with singularities, e.g., folding the previous planar interface along a straight line. In this 
way, we create a one dimensional straight singularity. Along the singularity, we can define 
a geometric solid angle. For example, let us consider a wedge of 90 degrees. In this case, 
the solid angle sustained by the surface at the singularity is 7r (Cl — 1/4), for the modes 
localized along the singularity, i.e. defined on the singularity only, the surface integral in Eq. 
(7) reduces to a line integral over the singularity. Any function defined along this straight 
singularity, only and anywhere else, will be solution of the line integral equation with the 
eigenvalue A = 0. Therefore using a Drude approximation of the dielectric constant of the 
metal and using Eq. (11), we obtain a plasmon wedge frequency u s = Ub/2. 

This derivation can be extended to a general wedge solid angle. Let it be a the interior 
angle defined by the intersection of two semi-infinite planes (see Fig. 2). In this case the 
solid angle sustained by the surface at the singularity is 2a (Cl = a/2ir). Using Eq. (11), 
we obtain the plasmon wedge frequency u s = u^a/27t. 




FIG. 2: Wedge geometry with an interior angle a. 
C. Modes localized at a sharp point 

Let us consider a point surface singularity, i.e., a metal vacuum interface with a sharp 
point. In this case, in Eq. (7), the integral over the singularity vanishes, giving also for the 
mode localized at the singularity, a A = mode as the only possible solution. 
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For the conical geometry, where the full metallic cone is defined by the apex angle 2a, 
illustrated in figure Fig. (3), we calculate using Eq. (11), the frequency of the mode 
localized at the apex. The geometric solid angle is given by fl s = 2n(l — cos(a)), giving for 
the frequency of the mode, called here a dot plasmon Ub mode: 



2 1 — cos a 2 
CO = z 0J h 



(12) 



It is interesting to notice that for small angle, this mode becomes linearly angle dependent, 
Ud ~ Uba/2. These results are relevant due to the important role of the tip geometry, for 
the tip-enhanced spectroscopies, when the tip can can be modelized by a metallic cone. 



a 



V 



FIG. 3: Conical geometry with an apex angle 2a. 

These results are generalizable to modes localized at corners. Let us consider the corner 
at the intersection of three semi-infinite planes normal to each other, where the metallic 
part is filling 1/8 of the space. We have at the intersection of the planes the solid angle of 
Q s = 7r/2. Using equation (11), we obtain the frequency of the mode localized at the corner, 

Ud = 27T 6 - (13) 

D. Complementarity rule 

One of the important properties of the plasmon modes demonstrated in smooth surfaces 
is the complementarity rule that relates the plasmon energies of a system and those of its 
complementary^. Here, we are going to demonstrate that it can be generalized to surfaces 
with singularities. The important point, as it is obvious from Eq. (7), is that the eigenvalues 
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Semi-infinite 
medium 


Wedge(90°) 


Wedge (a) 


Corner(90°) 


Dot plasmon(a) 
conical 


(uj/uj b ) 2 


1/2 


1/4 


a/2vr 


1/8 


(1 — cos a)/2 


tint/ text 


-1 


-3 


1 - 2vr/a 


-7 


(cos a + l)/(cosa — 1) 



TABLE I: Comparison for different geometries of surface singularities squared of the plasmon 
energy divided by the bulk plasmon energy within the Drude model (see Eq. 11). Values of the 
ratio tint I text-, between the interior and the exterior dielectric functions at the frequency of the 
surface plasmon for different geometries of surface singularities (see Eq. (10)). 

of a system X(S) and the ones of its complementary X(S) are the sameP 

X(S) = X(S). (14) 

The reason is that the two systems show the same surface, the only difference is that the 
values of the dielectric constant at both sides of the surface are interchanged. In other words, 
to treat the complementary system, we just need to interchange the values of e int and e ext 
in Eq. (10). For a Drude metal, one finally obtains 

^{S)+^(S)=ul (15) 

This rule is obviously applicable to any surface geometry, i.e., no matter it is a smooth 
surface or with a sharp singularity where a solid angle can be defined. 

E. Extension to general dielectric materials 

If the metal is not a Drude metal and/or the second medium is not vacuum, one can 
still get the plasmon frequencies using the more general equation(lO). In this case, one 
can get the plasmon frequencies using experimental data for the dielectric function of both 
systems^. 

In table I, we summarize for the different cases explained above, the results for the 
plasmon frequencies obtained for Drude metals, and the conditions to be fulfilled by the 
ratio of the dielectric function at the plasmon frequencies when the more general (10) has 
to be used. 
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IV. SUMMARY 



Herein, we have derived a general framework (7-9) and simple geometrical rules Eq. (10- 
11) which determines surface plasmon energies and plasmon modes, based on the shape of 
the metal. The results presented here provide simple and straightforward rules to design the 
energy of surface plasmons in several situations of experimental interest such as in plasmon 
wave guiding and in tip-enhanced spectroscopies. We have shown applications of the rule 
to obtain the surface plasmon energies in wedges, corners and conical tips. 

If the dispersion were included, one would be able to give valuable information for the 
plasmon wave guiding. Knowing that the framework of this theory allows to include itP^, it 
will be a further point of interest in the future. 



Appendix A: General derivation of a gradient formulation of the laplace integral 
equation 

Let us consider the partial differential equation 

L(w) = in B, (Al) 

where B is a region in 3D-space with surface S, and L is a linear elliptic differential operator 
operating on a sufficiently smooth function w of 3D-variables x. The Green's reciprocal 
identity for the operator L can be written: 

J B luL(v)-vL(u)] d B = J s (u a £-v^y S , (A2) 

in which u and v are arbitrary functions, L is the operator adjoint to L and n is the outward 
normal to S. Let us choose v as the fundamental solution G(x, x') = l/4-7r|x — x'| (in this 
case L = V 2 ) such that: 

L(G) = -<5(x - x') (A3) 

where 5 is the Dirac delta function. The integrands in Eq.(A2) can be expanded as an 
appropriate limit process^. We therefore consider a point x = x Q on the boundary S and let 
there be a sphere around this point separating a volume B p and surface S p from the rest of 
B. Let B denote what is left of region B with surface So comprised of S\ and S p according 
to the Fig. 4 




FIG. 4: Segment B p and the remainder Bq of the region B. 

Evaluating Eq.(A3) for B and So , which does not contain the point x , only the surface 
integral is left: 



(A4) 



The first term of Eq. (A4) has a singularity which is integrable, while the second term has to 
be treated with more care. Splitting the integral over So as S\ and S p and adding/ subtracting 
u(xq), and in the limit p — > 0, we get the main result iri^for xq on the boundary): 



c(x )tt(x, 



(A5) 



where the second integral is to be interpreted in Cauchy Principal Value sense. Obviously 
by setting u = lwe have: 



C(Xn 



47T 



(A6) 



where we have used that O(x ) is the solid angle pertaining to the point x . 
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